(12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(19) World Intellectual Property Organization 

International Bureau 

(43) International Publication Date 
15 February 2007 (15.02.2007) 




PCT 



(10) International Publication Number 

WO 2007/018810 Al 



(51) International Patent Classification: 

G06F 19/00 (2006.01) G01 V 3/08 (2006.01) 

(21) International Application Number: 

PCTYUS2006/025247 

(22) International Filing Date: 28 June 2006 (28.06.2006) 



(25) Filing Language: 

(26) Publication Language: 

(30) Priority Data: 

60/701,817 



English 
English 



22 July 2005 (22.07.2005) US 



(71) Applicant (for all designated States except US): EXXON- 
MOBIL UPSTREAM RESEARCH COMPANY; 

P.O.box 2189, Houston, TX 77252-2189 (US). 

(72) Inventor; and 

(75) Inventor/Applicant (for US only): LU, Xinyou [CN/US]; 
4028 Alexis Tate Circle, Missouri City, TX 77549 (US). 

(74) Agents: PLUMMER, J., Paul etal.; EXXONMOBIL UP- 
STREAM RESEARCH COMPANY, P.O.box 2189, Hous- 
ton, TX 77252-2189 (US). 

(81) Designated States ( unless otherwise indicated, for every 
kind of national protection available): AE, AG, AL, AM, 



AT, AU, AZ, BA, BB, BG, BR, BW, BY, BZ, CA, CH, CN, 
CO, CR, CU, CZ, DE, DK, DM, DZ, EC, EE, EG, ES, FI, 
GB, GD, GE, GH, GM, HN, HR, HU, ID, IL, IN, IS, JP, 
KE, KG, KM, KN, KP, KR, KZ, LA, LC, LK, LR, LS, LT, 
LU, LV, LY, MA, MD, MG, MK, MN, MW, MX, MZ, NA, 
NG, NI, NO, NZ, OM, PG, PH, PL, PT, RO, RS, RU, SC, 
SD, SE, SG, SK, SL, SM, SY, TJ, TM, TN, TR, TT, TZ, 
UA, UG, US, UZ, VC, VN, ZA, ZM, ZW. 

(84) Designated States ( unless otherwise indicated, for every 
kind of regional protection available): ARIPO (BW, GH, 
GM, KE, LS, MW, MZ, NA, SD, SL, SZ, TZ, UG, ZM, 
ZW), Eurasian (AM, AZ, BY, KG, KZ, MD, RU, TJ, TM), 
European (AT, BE, BG, CH, CY, CZ, DE, DK, EE, ES, FI, 
FR, GB, GR, HU, IE, IS, IT, LT, LU, LV, MC, NL, PL, PT, 
RO, SE, SI, SK, TR), OAPI (BF, BJ, CF, CG, CI, CM, GA, 
GN, GQ, GW ML, MR, NE, SN, TD, TG). 

Declarations under Rule 4.17: 

— as to applicant's entitlement to apply for and be granted a 
patent ( Rule 4. 1 7( ii) ) 

— as to the applicant s entitlement to claim the priority of the 
earlier application (Rule 4.17(iii)) 

Published: 

— with international search report 

For two-letter codes and other abbreviations, refer to the "Guid- 
ance Notes on Codes and Abbreviations" appearing at the begin- 
ning of each regular issue of the PCT Gazette. 



(54) Title: METHOD FOR DETERMINING RECEIVER ORIENTATIONS 



oo 




Create initial model for inversion 



1 



Invert data for receiver orientations 



o 



(57) Abstract: Method for completely specifying orientation of electromagnetic receivers dropped to the ocean bottom in an elec- 
tromagnetic survey. Survey data are selected, rejecting noisy data with long offsets and data where the receiver has saturated wit 
short offsets (61). A model is developed comprising three independent receiver orientation angles completely specifying the receiver 
orientation in three dimensions, and an earth resistivity model including a water layer and possibly an air layer (62). Maxwell's equa- 
tions, applied to the model and the selected data, are then inverted to determine the receiver orientations (63). 
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METHOD FOR DETERMINING RECEIVER ORIENTATIONS 

[0001] This application claims the benefit of U. S. Provisional Application No. 
60/701,817 filed on July 22, 2005. 

FIELD OF THE INVENTION 

[0002] This invention relates generally to the field of geophysical prospecting 
including reservoir delineation and, more particularly, to controlled-source 
electromagnetic surveying for geophysical applications. Specifically, the invention is 
a method for determining electromsagnetic receiver orientations. 

BACKGROUND OF THE INVENTION 

[0003] The marine controlled-source electromagnetic ("CSEM") exploration 
method uses man-made sources to generate electromagnetic (EM) waves and deploys 
receivers on the seafloor to record electromagnetic signals. The recorded 
electromagnetic signals are analyzed to infer subsea floor structures and/or determine 
nature of particular structures such as reservoirs. Figure 1 illustrates a typical 
deployment of CSEM equipment, with an horizontal electric dipole (HED) source 11 
towed by a vessel above the water bottom 12 on which receivers 13 are placed. This 
technology has been applied in hydrocarbon and mineral exploration, and also in other 
areas such as tectonic studies and environmental and geological engineering. 

[0004] At present,* receivers fall freely to the seafloor and therefore their 
orientations are unknown. Receiver orientations are required to determine the three- 
dimensional EM field vectors measured at receiver locations. The measured fields are 
then decomposed into components in preferred directions (for example, inline, 
crossline, and vertical) for analysis, inversion and interpretation. Effects on the 
decomposed components could be significant when the receiver cannot be oriented to 
those preferred directions because of inaccurate receiver orientations. Therefore the 
determination of receiver orientations could significantly affect data interpretation. 
The present invention provides a technique to determine receiver orientations. 
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[0005] In order to completely measure three-dimensional EM fields, receivers 
need be equipped with three mutually-perpendicular antennas for electric fields and 
three mutually-perpendicular magnetic sensors for magnetic fields. Three angles are 
necessary and sufficient to uniquely define the receiver orientations. These three 
angles establish the relationships between the measurement coordinates and receiver 
coordinates. A number of ways can be used to define the receiver orientations in the 
measurement coordinates. They are equivalent and can be converted one another. 
One way to define the receiver orientations is using azimuth and tilts for two 
horizontal channels (Figure 2). In Figure 2, (X, Y, Z) are assumed to be the 
measurement coordinates with X directed to the geodetic east, Y to the geodetic north, 
and Z upward. (X"\ Y">, Z'") are the receiver coordinates and the designed "east", 
"north" and vertical channels. (X\ Y\ Z') and (X", Y'\ Z") are auxiliary coordinates 
to help transform coordinates between (X, Y, Z) and (X"\ Y"\ Z'"). X 5 is the 
projection of X'" on the horizontal plane XY 5 while Y" is the projection of Y"\ 
With those setups, the receiver azimuth (a) is defined the angle between Y and Y 5 , the 
east channel tilt (J3) is the angle between X 5 and X" 5 , the north channel tilt (f) is the 
angle between Y' 5 and Y 5 ". 

[0006] At present, a number of methods are known for determining receiver 
orientations. One is to directly measure receiver orientations. Present direct 
measurement systems for marine CSEM receivers have reliability problems. Even 
with such a system available, an alternative method would be still necessary in some 
circumstances: for example, 1) when the direct measurement system is broken; 2) 
when the measurement system is not aligned with the directions of 
electrodes/magnetic sensors because of bending of long electric antennas on the 
seafloor. 

[0007] Another method is polarization analysis, disclosed in Behrens, J. P. (2005), 
"The Detection of Electrical Anisotropy in 35 Ma Pacific Lithosphere: Results from a 
Marine Controlled-Source Electromagnetic Survey and Implications for Hydration of 
the Upper Mantle," Ph.D. Thesis, University of California, San Diego (2005). Also 
see Constable and Cox, "Marine controlled source electromagnetic sounding 2: The 
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PEGASUS experiment/ 5 Journal of Geophysical Research 101, 5519-5530 (1996). 
The method is based on the fact that the EM field amplitude of the signal recorded by 
a receiver is maximized when the receiver antenna is in the direction of the transmitter 
(i.e. the major axis of the polarization ellipse) provided the transmitter is towed 
directly towards the receiver. Polarization analysis was the primary method used in 
early marine CSEM work to determine receiver azimuth in the subsequent data 
processing. The method requires at least one towline towed directly over each 
receiver. Receiver azimuth accuracy provided by this method is not very high. The 
average error in receiver azimuths is larger than 5 degrees for field data from a boat 
with a dynamic position system. It could be worse for a ship without a dynamic 
positioning system in rough weather conditions. 

[0008] Behrens also proposed to use coherency and correlation in natural EM 
signals recorded by different receivers to determine relative azimuth. This method 
was developed for receivers without a directly over-towing towline to complement the 
polarization analysis. The method determines the relative azimuth angle between two 
receivers. In order to find the receiver azimuth, the method requires the azimuth of the 
reference receiver be known. Success in using this method is dependent on whether 
high quality natural signals are recorded by both receivers. Accuracy by this method 
is normally lower than by polarization analysis. 

[0009] R. Mittet et al. used inversion to determine receiver azimuth in "E020: 
Inversion of SBL data acquired in shallow waters," EAGE 66 th Conference & 
Exhibition - Paris, France, June 7-10 (2004). This method overcomes limitations on 
both the polarization analysis and the method of using natural EM signals. All three 
of these methods, though widely used, address only the receiver azimuth, but do not 
disclose how to determine receiver orientations uniquely, i.e. both the azimuth and the 
tilts of the two horizontal channels. The reasons for neglecting the receiver's other 
two angles are at least three in number. (1) Data interpretation is mainly focused on 
(and data measurement may be limited to) the inline (meaning along the tow 
direction) electric component, which is normally not affected much by the tilts if the 
seafioor is not very steep. (2) The vertical electric component is either not measured 
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or is not fully utilized in data interpretation. (3) No reliable and accurate method is 
available to determine the receiver orientations. The two tilts are normally small ( < 
1 0 degrees because the seafloor is normally flat. The three reasons are obviously not 
completely independent of each other. 

[0010] In addition to the three reasons detailed above, the extent of possible 
impact on CSEM results of even small receiver tilt angles may not be appreciated. 
Effects of receiver orientations on the three electric components were simulated in the 
course of the present invention, and can be seen in Figs. 3-5. The source and receiver 
geometry used in the model calculations that generated these three drawings is taken 
from an actual field survey. The resistivity model is a layered earth model with water 
depth of 125 m. The towline direction is 265.57 degrees from the geodetic north, 
clockwise. In the modelings, the receiver (azimuth, a in Fig. 2) misalignment (da) 
with the towline ( the inline direction) is 1 5 degrees, the inline antenna tilt (fi) is up 5 
degrees, and the crossline tilt (y) is down 3 degrees. The modeling frequency is 0.25 
Hz. In each of these three drawings, the solid line represents an aligned and level 
receiver, the circles a level receiver with 8a = 15°, the + symbols an aligned receiver 
with tilts of p = 5° and y = -3°, and the broken line a misaligned and tilted receiver. 
Compared with the ideal situation (a level receiver aligned with the towline, i.e. all 
three angles [5a, J3 and y] zero), those figures show that while azimuth has much 
bigger effects on the two horizontal (inline and crossline) channels than do the tilt 
angles (especially on the cross component), the tilts have larger effects on the vertical 
component E z . The effects can be significant, for example, about one order in 
magnitude for the cross and vertical components of this example (Figs. 4 and 5). This 
example clearly demonstrates the importance of determining all three angles. 
Receiver azimuth alone cannot uniquely define the receiver orientations deployed on 
seafloor. 

[0011] In summary, there is a need for a technique to determine receiver 
orientations that can be used without any limitations on transmitter and receiver 
geometry. The present invention satisfies this need. 
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SUMMARY OF THE INVENTION 

[0012] In one embodiment, the invention is a computer-implemented method for 
determining three independent angles specifying orientation of electromagnetic 
receivers in a marine electromagnetic survey, comprising: (a) selecting survey data 
5 according to criteria including signal-to-noise ratio and degree of distortion; (b) 
creating a model representing the survey's source-receiver geometry and media for 
transmission of electric signal, said model comprising three receiver orientation angles, 
a resistivity model, and electromagnetic source (transmitter) parameters; and (c) 
solving Maxwell's electromagnetic field equations with the model and selected survey 
1 0 data as input information and said three orientation angles as unknowns, said solution 
being performed by iterative numerical inversion. 

[001 3] The invention is preferably practiced in the frequency domain in which case 
the survey data are transformed into the frequency domain by Fourier transform (or 
other method) before the selecting step above, or at least before the solving/inversion 
15 step. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0014] The present invention and its advantages will be better understood by 
referring to the following detailed description and the attached drawings in which: 

[0015] Fig. 1 illustrates a marine CSEM survey; 

20 [0016] Fig. 2 defines a set of three angles relating orientation of one coordinate 
system to another; 

[0017] Fig. 3 shows effects of receiver orientations on the amplitude of the inline 
electric field component; 

[0018] Fig. 4 shows effects of receiver orientations on the amplitude of the 
2 5 crossline electric field component; 
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[0019] Fig. 5 shows effects of receiver orientations on the amplitude of the 
vertical electric field component; and 

[0020] Fig. 6 is a flowchart of basic steps in one embodiment of the present 
inventive method. 

5 [0021] The invention will be described in connection with its preferred 
embodiments. However, to the extent that the following description is specific to a 
particular embodiment or a particular use of the invention, this is intended to be 
illustrative only, and is not to be construed as limiting the scope of the invention. On 
the contrary, it is intended to cover all alternatives, modifications and equivalents that 
10 may be included within the spirit and scope of the invention, as defined by the 
appended claims. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

[0022] This invention is a method for determining the orientation of an 
electromagnetic receiver in a marine CSEM survey by inversion of the 

15 electromagnetic field equations (Maxwell's equations). The three orientation angles 
are treated as inverted parameters, i.e., the unknowns to be solved for. The invention 
incudes three basic steps, summarized in the flowchart of Fig. 6: (i) prepare data for 
inversion (step 61); (ii) create an initial model (step 62); and (iii) invert data for 
receiver orientations (step 63). As used herein, inversion means inferring x from 

20 measured dataj/ = F(x, u, v, . . .), where finding x may be done by iterative trial and 
error using numerical methods such as finite difference techniques to find a solution 
that best fits the data by some selected criterion. 
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(i) Prepare data for inversion. 

[0023] Because of the well known skin depth effect, EM signals decay 
exponentially with distance from the source (or, transmitter) for a given frequency. 
The receiver cannot record high quality signals when the source is far away from the 
5 receiver because of ambient noises. When the source is too close to the receiver, the 
receiver is saturated because of the limited dynamic measurement range. So the 
measured signals are distorted. It is preferred for the invention that data are selected 
from such intermediate source-receiver offsets such that the source can generate 
signals strong enough at the receiver location to have good S/N (signal-to-noise ratio), 

1 0 but not too strong to saturate the receiver. In addition, accurate source and receiver 
geometry measurements are required for the selected data. The term "geometry" does 
not include angular orientation of the receivers, of course; as explained, this particular 
geometric feature cannot be measured with sufficient accuracy. The term does include 
orientations of transmitter (for example, azimuth and pitch for HED source) and 

1 5 coordinates of both the receiver and transmitter. 

[0024] Even within a selected offset range, data may not be ideal for inversion 
because of effects such as source instability, individual receiver electronic 
characteristic, temporally changing natural EM signals, and oceanic waves. The user 
of the present invention may wish to manually pick data to use, possibly with the help 
20 of interactive data display software, or according to experience. 

[0025] In CSEM surveying, both amplitude and phase are typically obtained for 
each EM field component that is measured. Either amplitude or phase data, or both, 
can be used for the inversion step of the present invention. For example, it could be 
that the phase data are assessed as having accuracy problem, in which case amplitude 
25 alone would be preferred for use in the inversion. In the most ideal situation, both 
amplitude and phase data of both the electric and magnetic fields, all six components, 
are included in the inversion. Some vertical component data are important to 
determine the tilts. Vertical component data are preferably about one third of the total 
data. Theoretically, at least 3 (independent) data points are needed to uniquely and 
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sufficiently determine the 3 angles of the receiver, where a value of Ex (either 
amplitude or phase) for one receiver/source position would constitute an example of a 
single data point. In practice, data for as many EM components as possible are 
preferably included because of noise and the different sensitivity of each component 
relative to each orientation angle. It is also preferable, but certainly not essential, to 
include as many frequencies and source-receiver combinations as possible. More data 
are more expensive to acquire, and require more computer time to process, but give 
more accurate results. 

[0026] CSEM survey data are measured in the time domain. The present 
inventive method is preferably performed in the frequency domain, in which 
embodiments the data must be transformed to the frequency domain by Fourier 
transformation or other methods. For example, the frequency content of the source 
waveform may be known, in which case amplitude and phase information for a 
specific frequency can be extracted from the measured data by data fitting techniques. 
All such methods shall be referred to herein as transforming the data to the frequency 
domain. In the frequency domain, the data become complex numbers. The present 
inventive method may be performed using only the real part of the selected data, or 
only the imaginary part, or both. Equivalently, as stated above, the invention may be 
performed with only amplitude data, or only phase data, or both. 
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(ii) Create an initial model 

[0027] An initial model is needed for inversion, which includes 3 receiver 
orientation angles and resistivity model. Electromagnetic source parameters such as 
source strength and frequency must also be included in the model as well as any 

5 needed (in the inversion step) receiver parameters such as receiver antenna length, and 
the field source and receiver geometry (the acquisition system must accurately record 
the geometry). This initial model should be created to be as realistic as possible. The 
inversion process (discussed below) must of course be performed by numerical 
methods, beginning with a first guess of the three receiver orientation angles. A good 

0 initial guess makes the inversion converge quickly and avoids the pitfall of a local 
minimum solution for the inversion. 

[0028] Receiver azimuth determined by other methods such as the polarization 
analysis can be used in the initial model. The seafloor slope is normally not steep, and 
therefore the angles of two horizontal channel tilts can be set to be zero in the initial 

5 model. The resistivity model can be a layered model consisting of air, seawater, and 
sedimentary seafloor. Seawater resistivity changes with depth and often is measured 
for each survey area. This measured seawater resistivity column should be used in the 
initial model if available. Otherwise, the seawater resistivity column can be estimated 
by exmpirical formula; see, for example, Chave et al, Electromagnetic Methods in 

o Applied Geophysics, M. Nambighian, Ed., Society of Exploration Geophysicists, Vol. 
2, 932 (1991). The sedimentary seafloor can simply be a half-space, or composed of a 
number of layers, or a more sophisticated model with inputs from other measurements 
such as seismic survey. 
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(iii) Invert data for receiver orientations 

[0029] The selected data are then inverted to determine the receiver orientations. 
Electromagnetic signals recorded by the receiver are related to transmitter and receiver 
geometry (both coordinates and orientations) as well as earth resistivity structures, 
5 This relationship (i.e. forward problem) can be written as d - F(m 9 r TX , o TX , r^ , ) , 
where d is a vector of measured data, m is a model of the earth's resistivity, 
( K rx >°Tx> r Rx>°Rx) are transmitter and receiver coordinates and orientations, and Fis 
known from Maxwell's equations for the EM fields and provides a means to compute 
d for any model m and transmitter and receiver geometry. It is well known that the 

10 earth's resistivity structures can be recovered from the CSEM measurements by using 
relationships provided by the above equation (Lu, et al. Geophys. J, Int. 138, 381- 
392(1999)). This process of solving the "forward" equation for certain variables that 
F depends on is called solving the inverse problem, or simply inversion. This 
invention uses the inverse process to recover receiver orientation. This inverse 

15 process can be simply written as follows: 

minimize object function = Wd - WGp + AR(m) , where W is a weighting matrix, 
G is a forward operator linearized from F 9 p is a vector of inverted parameters which 

includes both a model of earth resistivity and receiver orientations, R(m) is a 
regularization term to mitigate the nonuniqueness of inversion, and X is a 

20 regularization parameter. Acceptable answers may be obtained, however, with the 
regularization term set equal to zero, and using a least-squares iteration scheme. The 
inverse problem and its solutions have been studied extensively. See, for example, R. 
L. Parker, Geophysical Inverse Theory, Princeton University Press, Princeton, New 
Jersey (1994); W. Menke, Geophysical Data Analysis: Discrete Inverse Theory 7 

25 Academic Press, San Diego, California (1989); and A. Tarantola, Inverse Problem 
Theory, Methods for Data Fitting and Model Parameter Estimation, Elsevier, 
Amsterdam, The Netherlands (1987). In the course of developing the present 
inventive method, ID inversion software was licensed from University of Wisconsin- 
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Madison (Department of Civil and Environmental Engineering, Engineering School) 
and 3D inversion software was licensed from Saudia National Laboratory. 

[0030] The receiver orientation angles are needed in order to more accurately 
determine the subsurface resistivity structure. Accordingly, it may be difficult to 
5 arrive at a good guess for the resistivity model in the inversion for the orientation 
angles. In some embodiments of the invention, the prepared data are inverted for both 
the receiver orientation angles and the earth's resistivity model simultaneously. In 
other embodiments, the inversion is for the receiver orientation angles only, in which 
instances the results will depend on the accuracy of the assumed resistivity model. 
10 When the resistivity model is determined using the orientation angles found by one 
application of the present invention, the inversion for the orientation angles may be 
repeated, and then the resistivity model inversion can be performed a second time. 
This cycle may be iterated until desired stop criteria are obtained. The inversion 
calculations may be performed in ID, 2D, or 3D. 

15 [0031] Use of this method can determine receiver orientations. This invention 
uses inversion to determine all three angles which are necessary to define the receiver 
orientations, rather than just the receiver azimuth as proposed by Mittet, et al. 

Example 

[0032] A ID model with the same source and receiver geometry as the example of 
20 Figures 3-5 was used to generate a synthetic data set at frequencies of 0.125, 0.25, 0.5, 
1.25, and 2.0 Hz for a receiver with (a, fi 9 y) = (330.0, 2.0, 4.0). This data set was then 
used to test how well the present inventive method could determine the receiver 
orientations. In this example, the inversion process was set up for simultaneousely 
determining the receiver orientations and resistivity model. The initial resistivity 
25 model consisted of air, seawater and a uniform half-space for the sedimentary seafloor 
and the initial angles for receiver were (300.0, 0.0, 0.0). After applying the present 
inventive method, the resulting orientation angles were (a, /?, y) = (329.97, 1.95, 3.88). 
The recovered angles are very close to the angles used to generate the synthetic data, 
demonstrating the accuracy of the inventive method. 
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CLAIMS 

1. A computer-implemented method for determining three independent angles 
specifying orientation of electromagnetic receivers in a marine electromagnetic survey, 
comprising: 

5 (a) selecting survey data according to criteria including signal-to-noise ratio 

and degree of distortion; 

(b) creating a model representing the survey's source-receiver geometry and 
media for transmission of electric signal, said model comprising three receiver 
orientation angles, a resistivity model, and electromagnetic source (transmitter) and 

1 0 receiver parameters; and 

(c) solving Maxwell's electromagnetic field equations with the model and 
selected survey data as input information and said three orientation angles as 
unknowns, said solution being performed by iterative numerical inversion. 

2. The method of claim 1, wherein the survey data are transformed to the 
15 frequency domain at least before the solving by inversion step. 

3. The method of claim 2, wherein the transformation is by Fourier transform. 

4. The method of claim 1, wherein said source-receiver geometery includes 
angular orientations of the electromagnetic source (transmitter) and coordinates of both 
receiver and transmitter. 

20 5. The method of claim 1, further comprising: 

(d) revising the model to include the calculated receiver orientation angles 
and repeating the inversion step with resistivity as the unknown; and 

(e) repeating step (c) with the model's resistivity values updated to the 
resistivity values calculated in the preceding step. 
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6. The method of claim 1, wherein resistivity as a function of model position is 
treated as an additional unknown and is solved for in the inversion step. 

7. The method of claim 1, wherein the model is selected from a group consisting 
of a one-dimensional model, a two-dimensional model, and a three-dimensional model. 

5 8. The method of claim 1, wherein said selected survey data includes measured 
values of at least one inline electromagnetic field component, at least one crossline 
electromagnetic field component, and at least one vertical electromagnetic field 
component. 

9. The method of claim 2, wherein the selected survey data is selected from a 
l o group consisting of amplitude data, phase data, and both amplitude and phase data. 

10. The method of claim 1, wherein the inversion is performed by minimizing an 
object function. 

11. The method of claim 10, wherein said object function is 

object function - Wd - WGp 4- AR(m) 

15 where W is a weighting matrix, G is a forward operator linearized from Maxwell's 

equations, and p is a vector of inverted parameters which includes both a model of 
earth resistivity and receiver orientations, d is the selected data expressed as a vector, 
R(m) is a regularization term to mitigate nonuniqueness of inversion, and A, is a 
regularization parameter. 

20 12. The method of claim 1, wherein said selected survey data includes measured 
values of at least one inline electromagnetic field component and at least one vertical 
electromagnetic field component. 



25 



1 3 . The method of claim 1 , wherein survey data corresponding to source-receiver 
offsets greater than a pre-selected maximum value, and less than a pre-selected 
minimum value, are not selected for inversion, said maximum value being selected 
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based at least in part on signal-to-noise considerations and said minimum value being 
selected based at least in part on receiver saturation considerations. 

14. The method of claim 1 1, wherein said inversion is least squares inversion, and X 
= 0. 

15. The method of claim 1, wherein said selected data consist of at least three 
independent data points. 

16. The method of claim 2, wherein the selected survey data is selected from a 
group consisting of the real part of the data, the imaginary part of the data, and both 
real and imaginary data. 

17. The method of claim 1, wherein said resistivity model includes an earth layer 
and a water layer. 

18. The method of claim 17, wherein said resistivity model further includes an air 
layer. 

19. The method of claim 1, where said electromagnetic source and receiver 
parameters include source strength and frequency and receiver antenna length. 

20. A method for producing hydrocarbons from a subterranean region, comprising 

(a) obtaining electromagnetic field data from a controlled source 
electromagnetic survey of the subterranean region; 

(b) selecting survey data according to criteria including signal-to-noise ratio 
and degree of distortion; 

(c) obtaining a model representing the survey's source-receiver geometry 
and media for transmission of electric signal, said model comprising three receiver 
orientation angles, a resistivity model, and electromagnetic source (transmitter) and 
receiver parameters; 
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(d) obtaining solutions of Maxwell's electromagnetic field equations with 
the model and selected survey data as input information and said three orientation 
angles as unknowns, said solution being performed by iterative numerical inversion; 

(e) using the three orientation angles to resolve measured electromagnetic 
5 field data into inline, cross-line and vertical components, where line refers to a line 

containing successive source positions; and 



(f) producing hydrocarbons from locations in said subterranean region 
indicated by resistivity anomalies found in the resolved electromagnetic field data. 
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